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Abstract 

The nucleus of ^^^Thorium presents a unique isomer state of very low energy 
and long lifetime, current estimates are around 7.8 eV and seconds to hours respec- 
tively. This nuclear transitions therefore is a promising candidate for a novel type 
of frequency standard and severly groups worldwide have set out to investigate this 
system. Our aim is to construct a "solid state nuclear clock", i.e. a frequency stan- 
dard where Thorium ions are implanted into Calciumfluoride crystals transparent in 
vacuum ultraviolet range. As a first step towards an accurate determination of the 
exact energy and lifetime of this isomer state we perform low-resolution fluorescent 
spectroscopic measurements. 

Introduction 

Quantum frequency standards find a number of applications in fundamental science and 
applied technology. Nowadays the unit of time (SI second) is defined as 9,192,631,770 
periods of the radiation of the ground state hyperfine transition in ^'^^ Cesium [T]. Cesium 
fountain primary frequency standards have reached an accuracy level of 3 • 10~^^ [5]. 
At the same time, the progress in atom cooling and trapping together with sub-herz 
linewidth lasers allovifs to develop more accurate frequency standards based on optical 
transitions in neutral atoms and ions. For example, an accuracy level evaluation for a 
neutral ^^Sr optical lattice clock is 1.5 • 10~^^ [3], and the ^^Al+ quantum logic frequency 
standard demonstrated an inaccuracy level of order of 10"^'^ [3]. 

An interesting candidate for the role of an etalon transition in quantum frequency 
standards is the isomer transition in a nucleus of Thorium. This unique isotope has 
an extremely low-energy 229mfp]-^ isomeric state, expected at an optical wavelength of 
about 165 nm. The creation of a quantum frequency standards based on this transition 
will open new possibilities for experimental studies of basic fundamental laws of physics. 
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The question "are fundamental constants really constant" is among them. It is shown 
that the amplification factor between the drift of the nuclear transition frequency Suj/lo 
and the drift of the fine structure constant Sa/a has an absolute value between 10^ 
and 10^ [5]. Hence the Thorium quantum frequency standard will be several order of 
magnitude more sensitive to the possible drift of a than existing standards based on the 
transitions in electronic shell of atoms or ions. 

This paper is organized as follows. In section [1] we review current knowledge on 
the isomer transition energy and lifetime. In section [5] we briefly discuss two main 
experimental approaches to the direct observation of the isomer transition: trapped 
ions and Thorium-doped crystals. In section [3] we describe fluorescence spectroscopy 
of our Thorium-doped Calciumfiouride crystals and estimate the optical pumping rate. 
Section |4] gives an outlook on future steps. 

1 Isomer transition in the '^'^^Th nucleus 
1.1 Isomer state energy 

The isomer transition energy in ^^^Th nucleus is currently known only from undirect 
measurements. In 1976, Kroger and Reich [6] studied the 7-ray spectrum resulting from 
the decay of ^^'^ Uranium. They concluded that the ^^^Th nucleus has a = 3/2+ isomer 
state lying within 100 eV above the J'^ = 5/2+ ground state level. The development of 
high quality germanium detectors (resolution from 300 to 900 cV) allowed Helmcr and 
Reich to measure more precise 7-energies in 1989 ~ 1993 and to predict the energy of the 
nuclear transition to be £; = 3.5 ± 1.0 eV [7l[8]. 

This unnaturally low value triggered a multitude of investigations, both theoretically 
and experimentally, trying to determine the transition energy precisely, and to specify 
other properties of the J'^ = 3/2+ excited state of ^^^Th (such as lifetime and magnetic 
moment). However, direct searches for photon emission from the low- lying excited state 
have failed to observe a signal [SJ [TUJ [TTJ [H] . In 2005, Guimaracs-Filho and Helene [15] 
re-analysed the old experimental data and reported E ~ 5.5 ± 1.0 cV. 

New indirect measurements with the most advanced X-ray microcalorimeter (reso- 
lution from 26 to 30 eV) were performed by Beck ct.al in 2007 [M]. They published a 
new value for the transition energy E — 7.6 ± 0.5 eV for the isomer nuclear transition, 
shifting it into the vacuum ultraviolet domain. This shift probably explains the failure 
in observing the transition in previous experiments. In 2009 Beck et. al. reanalyzed 
their results taking into account the non-zero probability of the 42.43 keV — > 229mrp]_^ 
transition rate (estimated as 2%) and published a revised version E = 7.8 ± 0.5 eV [15j . 

The value E = 7.8 ± 0.5 eV is currently most accepted by the community but can 
not be considered definite untill the direct measurements are performed successfully. 
The dominant uncertainty in the prediction of Beck ct.al is connected with the value 
of the branching ratio b = 1/13 [14] from the 29.19 keV level to the ground state [16j . 
Estimation of this value vary in different works. Sakharov [16] mentions two alternative 
values for this branching ratio: 25 % in [17] which would result in £' = 9.3 ± 0.6 eV and 
51 % in [13] leading to E = 14.0 ± 1.0 eV. 
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1.2 Isomer state lifetime 

While the direct measurement of the 3/2+ — >■ 5/2+ ultraviolet transition in ^^^Th is still 
not performed, the isomer state lifetime r remains unknow. Most recent estimations of 
the half life of the isomer state in a bare nucleus based on theoretical calculations of the 
matrix element of the magnetic moment for the transition between ground and isomer 
levels were performed in jl8j . The theory was checked by comparison with experimental 
data for transitions at higher energies. They found the half-life of this transition T1/2 = 
(10.95h)/(0.025i?^), where E is given in cV. This corresponds T1/2 = 55 min for E = 7.8 
eV. 

For a Thorium atom or low charged ion the isomer lifetime can be significantly re- 
duced. In work of Tkalya [19] it is predicted that the probability of spontaneous magnetic 
dipole emission in a transparent non-magnetic dielectric medium with refractive index 
n is times higher than in vacuum. The reason for this factor is the renormalisation 
of the density of emitted photon states. For a CaF2 crystal, n = 1.55 at wavelength 
A = 159 nm (corresponding E = 7.8 eV), which gives a reduction factor of 3.7 for the 
lifetime. 

A more significant lifetime reduction can be connected with electronic bridge and 
bound internal conversion effects, i.e. isomer state decay can be amplified by interaction 
between the nucleus and the electron shell. There are a number of theoretical calculations 
of electronic bridge processes. In an isolated neutral Thorium atom the isomer state 
lifetime is about 10"^ s for E = 7.6 eV and about 4.5 min for E = 3.5 eV [5D]. In the 
Th^"*" ion for E — 7.6 cV the lifetime remains the same if the valence electron is in 
the ground state, and decreases 20 times if the valence electron is in the metastablc 7s 
state [21]. Similar calculations for the Th^+ ion which is most likely to occur in the 
solid-state approach have yet to be performed. In the Th+ ion the lifetime decreases 
10^ — 10^ times for E = 3.5 and E = 5.5 eV [12], there is not sufficient data about the 
Th+ ion spectrum to calculate the lifetime a.t E = 7.8 eV. 

Experimental attempts to measure the 229m'j']^ lifetime were performed in [23l 124) . 
These experiments are based on the theoretical prediction [26] that the a decay rate for 
229mrpj_^ is 2 to 4 times higher than for the ground state and that the energy spectrum of 
the a particles slightly differs for ground and isomer state. 

In [53] Thorium was electrodeposited onto a hollow cathode and then electronically 
excited by electric discharge. With some probability the excitation of the electronic shell 
was transfered to the nuclei [25] and therefore some amount of ^'^^'^Th was obtained. Af- 
ter the end of the discharge the a activity of the sample was monitored and a counts with 
different energies were written to different channels to record a spectrum. It was shown 
that a particles from Thorium are well energy separated from those of Thorium daugh- 
ter products. The excitation and detection procedure was repeated several times and a 
decaying component with half-life T1/2 = 2 ± 1 min was observed in counts for a parti- 
cles emitted from Thorium. In the same work also plots indicating a time dependence 
of count ratios for a particles from daughter products are presented [23] . These plots 
demonstrate some time dependence higher than 3(T variation (total number of counts 
is about 8000 per minute whereas the variation is about 300 per minute). This could 
indicate some technical artifact being responsible for the observed decay of the a count 
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rate from Thorium. 

In the work of Kikunaga et.al [24] a sample of 233u ^^^g chemically purified from 
Thorium, then it was left for Ih to allow the growth of ^^^Th in its ground and isomer 
state, then the Thorium was extracted. The a spectrum from this fresh Thorium sample 
was recorded during the first 6000 s and the following 6000 s. They compared the a counts 
in the region of 4915-4955 keV, where the a particles are produced in the decay of the 
229mfp]-^ isomer state with much higher probability than in the decay of the ground state. 
Slightly more a counts where detected for the first time interval and the authors derived 
< 2h with 3(7 confidence. This experiment was not suited to measure half-lifes on 
the order of minutes because of the long chemical preparation process. 

In general the analysis of these experiments is further complicated by the presence of 
various chemical compositions of Thorium with different decay rates depending on the 
chemical state (mainly ThO"*", Th and Th"*" in [23] and hydroxide and chloride complexes 
in [H]). So far, no clear evidence for a particles or photons originating from the 229mr[i]_^ 
isomer state has been presented. Therefore both the energy and the half-life of the 
229m.rYYy jgomer state are still not determined definitely. 



2 Approaches to direct measurements 

There are several groups worldwide working towards a direct measurements of the Tho- 
rium isomer transition. (T. Schumm, Institute of Atomic and Subatomic Physics TU 
Vienna, Austria; E. R. Hudson, University of California, USA E. Peik, PTE, Germany; 
M. Chapman and A. Kuzmich, Georgia Institute of Technology, USA). The expected 
excitation energy of 7.8 eV enforces work with Thorium in an ionized state to avoid over- 
whelming optical response by the electron shell which would mask the nuclear signal. 
Two main experimental directions are proposed. The first is the spectroscopy of Tho- 
rium ions in ion traps [571 UHl HH] ■ The second is the solid-state nuclear spectroscopy, 
where Thorium ions are implanted into a crystal (such as CaF2 or LiCaF) transparent 
for vacuum ultraviolet light [301 131] ■ 

The main advantage of solid-state nuclear spectroscopy is the possibility of simulta- 
neously excitating a macroscopic number (just a few orders of magnitude smaller than 
total number of atoms in a crystal) of Thorium nuclei whereas for ion trap spectroscopy 
this number is limited by 10^ 10^. This high number facilitates the first search for the 
transitions compared to ion traps. The second advantage is the fact that the experimen- 
tal setup does not require any bulky equipment for ion trap creation such as ultrahigh 
vacuum and laser cooling. However, the solid-state spectroscopy method is not free of 
additional difficulties connected with the hardly removable or irremovable systematic ef- 
fects in a crystal. Here, trapped ions present a very clean system with well-controllable 
parameters, probably ultimately more relevant for the realization of a high-performance 
laboratory frequency standard. 

In our study we plan to first perform low-resolution fluorescence spectroscopy using 
CaF2 as a host crystal to find the transition with 1 nm precision. CaF2 has a large (more 
than 10 eV) band gap. Th^+ ions can replace Ca^"*" ions in the crystal lattice (ionic radii 
are 94 pm for Th"*"*" and 114 pm for Ca^+), charge is compensated cither by a Ca vacancy 
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Figure 1: Scheme of fluorescence spectroscopy procedure. It lias three steps: a) irradia- 
tion of the sample by vacuum ultraviolet light; b) decay of own crystal fluorescence; c) 
count of the long-lived fluorescence 

or two F interstitials. 

CaF2 is not a good scintillator at the room temperature which is important as ^^^Th 
nuclei undergo a decay with a half-life of 7.8 • 10"^ years (a decay rate Fq = 2.8 • 10^^^ s). 
According to the total time of a photon scintillation burst following an a decay 
does not exceed 20 /is. The scintillation decay constant for events caused by 7 or X 
ray accompanying most of the a decays is 2 times longer therefore the total decay time 
tsc should not exceed 50 pjs,. ^^^Th has 8 short-living radioactive daughter elements. 
Therefore the number Usa of decays per second in the sample containing Nth Thorium 
nuclei with equilibrium number of daughter nuclei is about 

n,c = ^NThTc ^ 2.3 • IQ-^^Nth- (1) 

The scintillation will perturb the detection system during the part Usc ■ tgc of data aqui- 
sition time. The requirement rigc • i^c 1 gives the upper limits for the total number of 
Thorium nuclei in the sample: Nxh ^ 8.7 • lO^'' for tsc = 50 /iS. 

Another important question concerns the position of possibly emerging electronic 
states connected with Thorium impurities in the crystal. If there is a level with an 
energy near the isomer transition energy, it can increase the rate of isomer state decay 
via the electronic bridge process. This critical point is issue of ongoing theoretical and 
experimental work. 

3 Fluorescent spectroscopy 

Low-resolution fluorescence spectroscopy is the first step towards an accurate determi- 
nation of the exact energy and lifetime of the '^^^"^Th isomer state. Figure [1] shows a 
scheme of the planned spectroscopy experiment. We plan to irradiate the sample for a 
time texc with the probe radiation (laser or synchrotron) . We wait for a dead time to 
let the crystal fluorescence decay and then collect the fluorescence photons during time 
tdetect- Scanning the frequency of radiation and analysing the obtained photon counts, 
we plan to find the frequency of the isomer transition. 

A recent study by the Hudson group using the synchrotron radiation available at the 
Advanced Light Source showed that typical relaxation time of crystal fluorescence (not 
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connected to the isomer decay) for similar crystals (LiCAF, LiSAf NaYF etc.) doped with 
^■^^Th is less or about 1 s [33] . We plan to perform similar measurements for ^^^Th-doped 
CaF2 crystal using a high power laser system at the Max Born Institute in Berlin [M] ■ As 
long as the isomer state lifetime considerably exceeds the crystal fluorescence lifetime and 
the isomer transition lies within the band gap of CaF2 it should be possible to determine 
the isomer transition energy using the solid-state approach. 

The number of photons detected from a spontaneous nuclear isomer decay is: 

= (1 _ e-r*=..) X e-^^'' x (l - e"!^*---) , (2) 

where F is the total decay rate, W is the optical pumping rate, Fg is the radiative decay 
rate of the bare nucleus, n is the refraction index, fi is the solid angle covered by the 
detector, k is the detector efficiency. 

We now estimate the optical pumping rate W. For the low-resolution spectroscopy 
stage, the spectrum of the radiation source will be broad and smooth, therefore W 
depends only on the spectral power density dl/du of the radiation at the transition 
frequency uj = ujeg- Let us consider for the sake of simplicity the Lorentzian spectral 
density of intensity: 

rf/ ^ dh 1 

duj dcj 1 + 4(u; - cjeg) ' 

where Tl is the spectral linewidth of the source, ^ is the spectral intensity density on 
resonance. 

Th^+ ion has a Rn-like electronic shell (noble gas) with zero angular momentum and 
zero spin. In the following we neglect the structure of the electronic shell as well as 
higher levels of the Th nucleus. Ground and isomer state are degenerate with respect 
to projections m and m' of angular momenta J and J' respectively (here we neglect the 
interaction with ramdomly oriented spins of surrounding F nuclei) . The magnetic field 
B of the spectroscopy radiation can be represented as B = Bo cos(wt -I- f{t)), where ip{t) 
describes the random phase variations responsible for non-zero spectral linewidth F^. 
Note that for the low- resolution spectroscopy F^ just a few orders of magnitude smaller 
than UJ. 

The Hamiltonian H in resonance approximation can be represented as 

H = Ho + hJ2Veg [|g)(e|e*("*+'^W) + |e) (g|e-^('"*+'^(*»j + H' . (4) 

Here Veg = BQqfi1g/{2h), B^q is the gth covariant component of Bq, /i^^ is the g-th 
contravariant component of the matrix element fl^q of the nucleus magnetic momentum 
operator in cyclic coordinates, q = m' — m. Hq is the bare Hamiltonian of the Thorium 
nucleus, H' describing the interaction of the nucleus with all other fields except the spec- 
troscopy radiation. We can choose a quantization axis, for example, along the radiation 
propagation direction, and classify excited states e and ground states g by angular mo- 
mentum projection onto this axis. The set of equations for the density matrix pij in the 
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rotating frame with relaxation terms reads as: 

Pee = X! ^Pa" ~ P^s^ ^ ^P"" + X! -T P'''"' ' 

a e' 

pee' = -i ^ [VegPge' - Ve'gPeg] - (T + r']_)pee' , 6 7^ e'; (6) 
9 

Peg = - ^')Peg - i ^ VegPg'g + i ^ VJj'gPee' - F'; (7) 
S' e' 

- XI ^^=3 [^-=9 " ^S"^] + XI e^egPee + ^ Ps'ff' ; (8) 



Pas 



a 



Pgg' = -i^iVegPeg' -Veg'Pge]~T'lPgg,, 9 ^ g' ■ (9) 

e 

Here F is the decay rate of the isomer state, and F||"' are decoherence and mixing 
rates between magnetic states of the same energy level, F' is the relaxation rate of optical 
coherence p^g- We can estimate F^^ and F'|-' via the spin-spin relaxation time r ~ 20/xs 

for CaF2 [35]. Therefore, F||-' « T^^^ ?» l/r « 5-10^ s^^. The optical coherence relaxation 
rate F' can be represented as F' = Fi/2 + F/2 + F^f, [35], where Tab is the relaxation 
rate due to adiabatic perturbation. For any reasonable low-resolution radiation source 
we can take 

Kg < F <F||^ «F^^' < F' = Fl/2. (10) 
From (5-10) and the normalization condition pa = 1 we find: 

5 ' 



1 IV 

Pe.c = W-Tp,^,- M/=-X^7^- (12) 

where pexc = ("ee IS the total excited state population. 

Now we should calculate X^egKg- For the sake of simphcity we suppose that the 
radiation has a certain polarization along q'th cyclic unit vector, so B^q = \BQ\5qqi. 
In dipole approximation, F — X]g l/^egP- According the Wigner-Eckart theorem, 

Meg = \\^^\\^Jm^q^ whcrc \\p\\ \s a reduced dipole moment, Cj™iq is a Clebsch-Gordan 

, , 2 

coefficient, c is the light velocity in vacuum. Using „ C"' ™ 



^ Jrnlq 



1, we find 



Hl = \/^- (13) 



On the other hand. 



_ Bi\\p\\' ^ I ^_ i?g||/i|P(2J^ + l) . . 

^eg I^Jmlg' 3 ' ^-^ 

e,g m,m' 
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Tabic 1: Estimation of optical pumping rate W for diffcrnt sources 



Source 


A, nm 


dl/dX, W/cm2-nm 


W, s-i 


McPherson 632 thermal lamp 


113-200 


2.5 • lO-'^ (165 nm) 


6 • 10-16 


Q-switched laser (NT342B-SH) 


210-420 


4 (225 nm) 


10-9 


Synchrotron (ALS, USA) 


40-250 


32 (165 nm) 


3.25-10-9 



The amplitude |i?o| of magnetic field is easy to express via the total radiation intensity 
/: |i?o| = SttJ/c. Integrating (3) over w, we obtain / = '^^'^^ therefore 

From (12)-(15) we finally obtain the expression for the optical pumping rate W: 

_ 2 c^Ttt^ dip _ XT dip 

3 hio^ duj 2ATT^c'^h dX' ^ ' 

Here we took into account the values of total angular momentum J — 5/2 and J' = 3/2 
in the ground and excited state respectively. Typical values of W are estimated in Table 
1. Wc see that if the sample contains Nth = 10^'^ Thorium nuclei, and if we use a 
powerful excitation source (e.g. ALS synchrotron), we can excite up to 3.25 ■ 10'^ nuclei 
per second. 



4 Future plans 

We are currently in preparation for the low-resolution spectroscopy stage. Our immediate 
plan is to develop and improve noise reduction techniques and to measure the fluorescence 
rate of a ^^^Th-doped CaF2 crystal after irradiation with a high-power laser field to 
accurately determine the crystal fluorescence time r^j. In parallel we perform absorption 
measurments on ^^^Th-doped to determine whether addtional electronic states within the 
band gap are created by the doping process. Theoretical simulations of the microscopic 
properties of the doping complex are under way. These measurements will help us to 
determine the optimal measurement strategy concerning the timings texc, Td, and t detect 
as well as the required excitation source and the optimal doping concentration. 

In a second step we plan to grown the ^^^xh-doped CaF2 crystal and perform the 
low-resolution spectroscopy measurements as discussed above. 

Once the isomer energy will be found, we will proceed to the next step, the high- 
resolution spectroscopy. We plan to use high-harmonic generation (5th harmonic of 
800 nm) of a femtosecond frequency comb system to generate a narrow-linewidth ex- 
citation source with absolute frequency callibration. At this stage, we will be able to 
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quantify line shifts and broadenings due to the crystal environment and validate the 
general experimental approach. 

Finally we plan to actively stabilize the frequency comb to the nuclear transition and 
thereby realize to a solid-state nuclear frequency standard. 
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